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Abstract: We predict the amplitude of the B — > K£ + £~ decay in the region of the 
dilepton invariant mass squared < q 2 < m 2 ,^, that is, at large hadronic recoil. The 
B —7- K form factors entering the factorizable part of the decay amplitude are obtained from 
QCD light-cone sum rules. The nonlocal effects, generated by the four-quark and penguin 
operators combined with the electromagnetic interaction, are calculated at q 2 < 0, far below 
the hadronic thresholds. For hard-gluon contributions we employ the QCD factorization 
approach. The soft-gluon nonfactorizable contributions are estimated from QCD light-cone 
sum rules. The result of the calculation is matched to the hadronic dispersion relation in 
the variable q 2 , which is then continued to the kinematical region of the decay. The overall 
effect of nonlocal contributions in B — > K£ + £~ at large hadronic recoil is moderate. The 
main uncertainty of the predicted B — > K£ + £~ partial width is caused by the B — )• K form 
factors. Furthermore, the isospin asymmetry in this decay is expected to be very small. 
We investigate the deviation of the observables from the Standard Model predictions by 
introducing a generic new physics contribution to the effective Hamiltonian. 
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1. Introduction 

After a very successful start, the LHCb experiment has already provided new measurements 
of exclusive flavour-changing neutral current (FCNC) decays of B mesons [|}-[||], continuing 
the studies of these decays carried out at l?-factories [||, |6| and Tevatron 0] . The favourite 
decay channel is B — > K*£ + £~ where a rich kinematics allows one to measure several 
nontrivial observables sensitive to the underlying FCNC b —> s£ + £~ transitions. So far, 
the measurements of these decays did not reveal any significant deviation from Standard 
Model (SM). In particular, the zero crossing point of the forward-backward asymmetry in 
B —7- K*£ + £~ predicted in SM has recently been observed within the expected interval 
of the variable q 2 , the invariant mass squared of the lepton pair. Meanwhile, the measured 
isospin asymmetry in B — > K^£ + £~, allows for the values larger than expected in SM, 
hence, this observable deserves a further careful study. 

The theory of rare semileptonic decays, such as B — > K^£ + £~~ , suffers from hadronic 
uncertainties. Apart from the heavy-to-light form factors, determining the contributions 
of the leading FCNC operators in the effective Hamiltonian, the decay amplitudes receive 
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contributions of the current-current and penguin operators combined with the lepton-pair 
emission via electromagnetic (e.m.) interaction. A common feature of all these effects is 
that they are nonlocal because the quark-flavour transition is separated from the virtual- 
photon emission, and the characteristic distances of this separation are not necessarily 
small. The most important nonlocal contribution is generated by the current-current oper- 
ators with c quarks, forming "charm-loops" after photon emission. The charm loops turn 
into intermediate charmonium states at q 2 > m 2 ,^. The nonlocal effects in B — > K^£ + £~ 
have been estimated in the literature, within different models and approximations. The 
most complete analysis at large hadronic recoil (low q 2 ) was done in the QCD-factorization 
framework in H. In that approach, the quark-gluon diagrams of various factorizable and 
nonfactorizable nonlocal contributions were calculated at timelike q 2 , below charmonium 
thresholds. 

In our previous paper Q the soft-gluon emission from the charm loop in B — > K^£ + £~ 
- an effect not accessible in QCD factorization approach- was estimated. To that end, the 
light-cone operator product expansion (OPE) at q 2 <C Am 2 was employed and the soft- 
gluon emission was effectively resummed in a nonlocal quark-antiquark-gluon operator. 
The resulting hadronic matrix elements were then estimated using QCD light-cone sum 
rules (LCSR), i.e., the same approach which is used to calculate the B — > K^*' form 
factors. Being suppressed by the powers of A-q CD / '(Am 2 — q 2 ), the soft-gluon contribution 
to the decay amplitude grows at q 2 approaching the cc threshold, To avoid the divergence 
at q 2 ~ Am 2 , the hadronic dispersion relation in the q 2 channel was employed in |J and 
matched to the calculation result at q 2 <C Am 2 . 

The aim of this paper is to perform a complete account of hadronic nonlocal effects 
for one particular channel B — > K£ + £~. We will apply the same approach as in ||, this 
time including the complete effective Hamiltonian. In this case, light-quark loops also 
contribute to the amplitude, therefore we calculate the nonlocal effects in QCD in the 
spacelike region q 2 < 0, sufficiently far from all quark-antiquark thresholds. Soft-gluon 
emission evaluated in || is taken into account for quark loops with different flavours. In 
addition, we calculate a new effect of soft-gluon emission from the gluon-penguin operator 
applying the LCSR approach. We also take into account the hard-gluon NLO contributions. 
In the correlation functions used to derive LCSR's, the hard-gluon exchanges generate 
multiloop and multiscale diagrams which demand dedicated calculational efforts that are far 
beyond our scope. Instead, we approximate the hard-gluon nonlocal effects, employing the 
QCD factorization approach || at q 2 < 0. Following Q, the physical region in B — > K£ + £~ 
is accessed via hadronic dispersion relations in the variable q 2 . In these relations, in 
addition to the charmonium states, the light vector mesons contribute. Finally, we predict 
the observables for B — > K£ + £~ including the differential width and isospin asymmetry. 
Our results are applicable at q 2 < fn 2 j^ where they are compared with the available 
experimental data. 

The B —?• K£ + £~ decay channel chosen here has less observables and a smaller branch- 
ing fraction than B — > K*£ + £~. In turn, the relative simplicity makes the kaon mode a 
more convenient study object for hadronic effects. First of all, the current accuracy of the 
B —7- K form factors is better than for B — > K* form factors. Note that the latter form 
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factors are available only in the quenched approximation of the lattice QCD. The B — >• K 
form factors used here are calculated from LCSR with kaon distribution amplitudes (DA's). 
These sum rules are free from the "systematic" uncertainty of LCSR's for B — > K* form 
factors caused by neglecting the K* — > Kn width. The alternative LCSR's with S-meson 
DA's implicitly overcome this problem if one relies on the duality approximation in the K* 
channel. However, the i?-meson DA's still suffers from large uncertainties of their param- 
eters, such as the inverse moment. Since the form factors not only determine the leading 
operator contributions but also enter the factorizable nonlocal effects, we expect that the 
hadronic input in B — > K£ + £~ is currently under a better control than in B — > K*£ + £~ . 
This circumstance strengthens the case for using the kaon channel as a perspective tool to 



trace generic new physics (see e.g., []10|, pT[|). 

The plan of the paper is as follows. In Section 2 we define the nonlocal hadronic 
matrix elements entering the B — > K£ + £~ amplitude. The nonlocal effects are listed in 
terms of quark-gluon diagrams. In Section 3 we present the corresponding hadronic matrix 
elements in the spacelike region obtained by combining LCSR's with the QCD factorization 
approach. In Section 4 the LCSR for the soft contribution of the gluon-penguin operator 
is derived. Furthermore, in Section 5 we specify the use of hadronic dispersion relations 
to access the physical region. In Section 6 the relevant numerical analysis is performed. 
The resulting predictions for the B — > K£ + £~ differential width and isospin asymmetry 
are presented in Section 7 where we also discuss the influence of generic new physics 
on these observables. Section 8 is reserved for the concluding discussion and outlook. 
In Appendix A, the operators of the effective Hamiltonian and their Wilson coefficients 
are specified, Appendix B contains the definitions of S-meson DA's and Appendix C the 
expressions related to the LCSR derivation. 



2. Hadronic effects in the decay amplitude 

The B K£ + £~ decay amplitude in SM is given by the matrix element of the b —> s£ + £ 



effective Hamiltonian [12, 13 1: 



H eff 



4G F 



^VaK'EQWOiM . (2.1) 

* i=l 



The effective local operators Oi and the numerical values of their Wilson coefficients d are 
presented in App. A. In this paper we generally neglect the CKM-suppressed contributions 
proportional to V u bV* s , retaining them only in the calculation of the isospin asymmetry 
where their numerical impact is non-negligible. 

In the decay amplitude we separate the dominant factorizable contributions of the 
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(a) 

Figure 1: FCNC contributions to B 
denoted as black squares. 



(b) 

K£ + £~ due to the effective operators 09,10 (&) o-nd O7 (b) 



operators 07,9,10 (Fig. U): 



A(B -»■ K£ + t 



-(K(p)£+r \H eff \B(p + q)) 



V2 vr 



rriB + m K 



eff 

7. 



(2.2) 



where f BK (q 2 ) is the usual B — > if vector (tensor) form factor and, according to [13], 



a 



*ff - Cr - lC 5 



Ca. The nonlocal hadronic matrix element 



Hf^=i [ d A xe^^K(p)\T{jf l m (x),[c 1 O ( C\0) + C 2 O { 2 c) (0) 



+ Yl C k O k (0)]}\B(p + q)} = [(p. q )q^-qV]'H (BK Hq 2 ) 
fc=3-6,8 



(2.3) 



contains the contributions of all remaining operators in ( |2.1| ) combined with the quark 
electromagnetic (e.m.) current = Ylq=u d s cbQqQ'~1v c L- Our definition of T~L^ BK ^ differs 
from the one used in Q by including the quark-charge factors Q q . Furthermore, due to 

yeff 



renormalization effects, one has to replace above Cgg— >Cg JJ = Csg+i^C^—C^) /3. After 
substituting ( pj| ) in ( |2.2[ ), due to conservation of the leptonic current, only the structure 
q 2 p^ remains in (^3]), canceling the photon propagator 1/q 2 . Hence, opposite to the case of 
B —7- K*£ + £~ , there is no kinematic enhancement of the B — > K£ + £~ amplitude at low q 2 . 
Dividing the invariant amplitude T~L^ BK \q 2 ) by the form factor fg K (q 2 ), it is convenient to 
represent the nonlocal hadronic effect in a form of a (process- and g 2 -dependent) correction 
to the coefficient Cg, 



AC, 



(BK) 



l^ 2 U( BK )(q 2 ) 



(2.4) 



The contributions to the nonlocal amplitude fl2.3|) are usually represented in a form 
of quark-gluon diagrams shown one by one in Figs.|2| - ||. The dominant contribution 
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(a) (b) (c) 

Figure 2: Factorizable quark-loop contributions to B —¥ K£ + £~ due to four-quark effective op- 
erators 0\ 2 and 03_6- Crossed circles denote the possible points of the virtual photon emission. 
Diagrams similar to (c) with the gluon line attached to the s-quark line are not shown. 




(a) 



(b) 



Figure 3: Nonfactorizable quark-loop contributions to B — > K£ + £ : (a) with soft gluon (denoted 
by crossed line) and (b) with hard-gluon. 



is generated by the operators 0\ 2 with large Wilson coefficients and corresponds to the 
diagrams shown in Figs. [2], ||, where q = c (the "charm-loop" effect). This contribution 
includes intermediate vector charmonium states in the upper part of the decay kinematical 
region, ffiju < q 2 < (tub —rux) 2 - Measuring B — > K£ + £~~ in the vicinity of q 2 = (ifi = 
J/tp, ip(2S), ...), one practically observes a combination of the weak nonleptonic {B — > tpK) 
and leptonic (ip — > l + £~) decays, whereas the genuine FCNC process b — > s£ + £~, driven 
by the effective operators 07,9,10, turns into a tiny background. Note that the interval 
Am 2 D < q 2 < (ma — mx) 2 also contains broad charmonium resonances. To avoid the 
charmonium background, all experimental measurements of B — > K£ + £~ implement a 
subtraction of the two g 2 -bins around ra^ and in^^s) • 

In this paper we concentrate on the region of small and intermediate lepton-pair masses, 
q 2 < vt?ji^. In this region the charm-loop effect is present in the form of a virtual cc 
fluctuation. The latter is usually approximated by c-quark loop diagrams, starting from 
the leading-order (LO) simple loop in Fig. Q(a) and including the gluon corrections. A 
characteristic feature of this perturbative approximation at timelike q 2 is a "kink" in the 
predicted differential width, due to the onset of the imaginary part of the c-quark loop 
at the threshold q 2 = Am 2 < rn 2 ^. To avoid this unphysical effect, in what follows we 
employ the hadronic dispersion relation in the q 2 variable, following ||. In this paper we 
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(a) 



(b) 



Figure 4: Contributions of the Og g operator to B — > K£ + £ : (a) factorizable, (b) and (c) nonfac- 
torizable with soft and hard gluon, respectively. 



also include the contributions to 7i^ BK \q 2 ) stemming from the operators with light quarks. 
Hence, to stay away from quark-antiquark thresholds, all quark-gluon diagrams shown in 
Figs. ^-[| have to be calculated at q 2 < with sufficiently large \q 2 \. 




Figure 5: Weak annihilation contribution to B — > K£ + £ amplitude: (a) in LO and (b) one of the 
NLO hard- gluon exchange diagrams. 



Before considering separate contributions to the hadronic matrix element T~L^ BK \q 2 ), 
the following comment is in order. A close inspection of the quark-loop diagrams at NLO 
shows that also at q 2 < 0, below all hadronic thresholds in the e.m. current channel, the 
diagrams develop an imaginary part. This is a clear signal that the intermediate hadronic 
states in the amplitude H^ BK \q 2 < 0) go on shell, so that the quark-gluon diagrams 
provide only a local duality approximation for this amplitude. For example, let us consider 
a diagram similar to the one shown in Fig. ||(c), but with the photon emitted from the 
■s-quark. This diagram and one of its hadronic counterparts are shown in Fig. ||. To 
describe this hadronic transition in more detail we choose q = c, q' = d for definiteness. In 
this case 5°-meson decays to an on-shell hadronic state with the quark content (cs)(cd) 
(e.g., B° — > Dg D + ) which then converts via strong and e.m. interaction into the final 

meson and a lepton pair (D s D + — > K°£ + £ ). This hadronic mechanism involving 
multiple intermediate states with the same quantum numbers contributes to the strong 
final-state interaction in B — > K£ + £~ . Approximating the sum over intermediate hadronic 
states with a quark-level diagram is justified by the large mass of the initial S-state and 
the large recoil of the final hadron. This approximation is similar to the one used in the 
QCD factorization approach to exclusive nonleptonic B-decays [14]. 
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Figure 6: The NLO quark-loop diagram of B Kll decay (left) and one of its hadronic counter- 
parts (right) containing the intermediate pair of on-shell hadrons with flavour content (sq)(qq f ). 



Let us emphasize another important point. The imaginary part of the amplitude 
n^ BK) (q 2 ) at spacelike q 2 is generated by the discontinuities in the variable (p + q) 2 . In 
B — > K£ + £~ with an on-shell initial B meson, this variable is fixed at (p + q) 2 = m 2 B . 
To illustrate a generic two- variable kinematics, one has to consider a correlation function 
where, instead of the 5-meson state in (^), a vacuum state is taken and a quark current 
<?'75 b with B q i -meson quantum numbers is added to the time-ordered operator product, 
with the momentum squared (p + q) 2 . The variable q 2 is kept fixed and spacelike. Still, 
the correlation function has discontinuities in the variable (p + q) 2 located below the B- 
meson pole and generated by intermediate (sq)(qq') states. Note that the presence of these 
"parasitic" intermediate states in the 5-meson channel of the correlation function is caused 
by the absence of external 4- momentum in the effective operator vertex in (2.3). Returning 
to the calculation of the hadronic amplitude fi^ BK \q 2 ), we conclude that using OPE near 
the light-cone, x 2 ~ in ( |2.3] ) and restricting the calculation to q 2 < 0, one cannot avoid 
using quark-gluon diagrams at a timelike (albeit large) value of the second kinematical 
variable (p + q) 2 . Hence, strictly speaking, the calculation procedure used here relies not 
entirely on OPE but also on the local quark-hadron duality. 



3. Hadronic matrix elements in the spacelike region 



In this section we present separate contributions to the hadronic amplitude 7iBK(q 2 < 0) 
defined in (2.3), including all effective operators entering H e ff and taking into account 
different possible quark topologies. 



3.1 Factorizable loop 

In LO, the contribution of the four-quark operators to B — > K£ + £~ is shown in Fig. 0(a). 

(c) 

In addition to the dominant c-quark loop from the current-current operators 0\ l± it also 
contains quark loops with various flavours, originating from the quark-penguin operators 
O3-Q. In the adopted convention for the effective operators the sum over all LO loops can 
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be written as 



nj(BK) . 2s 
H fact,LOW ) 




+ C 2 g(m 2 c ,q 



+ 



c i + c t 



1 

8vr2 | - V 3 
Q s g(m 2 s ,q 2 ) + Q b g(ml,q 2 ) 



^2 Qq9{m 2 q ,q 2 



(3.1) 



q=u,d,s,c,b 



where 



g(m 2 q 2 ) 




(3.2) 



is the well-known loop function, n is the renormalization scale and g(m 2 , q 2 ) = g(m 2 ,q 2 ) + l. 
The definition of g(m 2 , q 2 ) used in ||] differs from the above by the factor 4/9. Since here 
we employ the loop diagrams at q 2 < 0, that is below all quark-antiquark thresholds, only 
the real part of g(m 2 ,q 2 ) is needed. Hereafter we also neglect the mass of the u and d 
quarks, so that for them the loop function is g{0,q 2 ) = 2/3 — ln(— q 2 //j?). 

3.2 Factorizable NLO contributions 

The NLO corrections to the quark loops include the diagrams shown in Figs. Q(b,c). In the 
same order of perturbative expansion the gluon-penguin operator enters with the diagrams 
in Fig. |](a). The characteristic feature of all these diagrams H 1 is that the correspond- 
ing hadronic matrix elements are factorizable in the same form as the LO contributions 
considered in the previous section. Following ||, we extract the hard-scattering coefficient 
functions multiplying the B — > K form factor from the NLO diagrams for the b — > s£ + £~ 
decay rates calculated in [[0], 17]. 

Summing over all contributing operators, the factorizable NLO contributions to the 
hadronic amplitude can be written as 



rt fact,NLO\ ( i ' 



a s mb 
327T 2 mg 



+ 



m B 
2mb 



Cl F?\q 2 ) + 2C 2 (if V) + ^ 2 ( V)) + C^'jtfV)] }/. 



<bk(Q 2 )- (3-3) 



The definitions and nomenclature of the indices of the functions ^128 are ^ ne same as m 
1 1 6 , |l^], where F^ 7 ^ are expressed as a double expansion in s = q 2 /m 2 and m 2 = m 2 /m 2 . 
It has been shown in |l6|, 17] that keeping the terms up to the third power of s and m 2 
provides a sufficient numerical accuracy in the region 0.05 < s < 0.25. Here we use these 
expansions at q 2 < 0, restricting ourselves to |s| < 0.25. For F^ 7 ' 9 ^ we use the expressions 
derived in l|. Note that the NLO diagrams generate an imaginary part in the coefficient 



1 For an earlier estimate of this effect within quark model see, e.g. (l£ 
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functions entering ( |3.3[) , which, as discussed in the previous section, has to be interpreted 
as a quark-hadron duality counterpart of the strong rescattering in the hadronic amplitude. 

A dedicated LCSR calculation of the hadronic matrix element ( |3,3| ) at q 2 < is also 
possible. One has to introduce a correlation function where the B meson is represented with 
its DA's and K meson is interpolated with a quark current. However, in the case of four- 
quark operators, the complexity of the two- loop diagrams with several momentum/mass 
scales emerging in the OPE for the correlation function makes the computation of hard- 
gluon contributions from LCSR's a very difficult task. On the other hand, the diagrams 
for the correlation function with the gluon-penguin operator have only one loop and can in 
principle be calculated in future, to be compared with the QCD factorization results which 
are used here. 

3.3 Nonfactorizable soft-gluon contributions 

The contributions to the hadronic amplitude T-Lbk(,Q 2 ) that are beyond QCD factorization, 
include also the soft-gluon emission from the quark loop, as shown in Fig. ||(a). In [6j the 
dominant part of this effect, caused by the c-quark loop was estimated using the light- 
cone OPE at q 2 <C 4m 2 , where this power correction is suppressed by the inverse power 
of (Am 2 — q 2 ) with respect to the LO loop. Here we also include the contributions of the 
soft-gluon emission from the quark loops generated by the quark-penguin operators. In 
presence of the light-quark loops, we shift the calculation to q 2 < 0, so that \q 2 \ 3> Aq CD . 

The B — > K hadronic matrix element for the soft-gluon emission from the loop with 
a generic quark flavour q can be easily derived from the results of ||. At q 2 <C 4m 2 , the 
soft-gluon emission effect is reduced to a hadronic matrix element 

(K( P )\6^m q ,q)\B(p + q)} = [(p • q)q„ - q%] A(m 2 q , q 2 ) , (3.4) 

where the operator 0^(m q ,q) is equal to the nonlocal effective operator 0^(q) defined in 
(3.14) of H with the replacement m c — > m q in the coefficient function I MpQ/ g(g,u;) given 
in (3.15). Collecting the contributions from all four-quark operators, we obtain for the 
soft-gluon emission effect 

^ioR^ 2 ) = 2{Q c (Ci + C 4 - C 6 )A(m 2 c , q 2 ) + (C 4 - C 6 )(Q U + Q d )A(Q, q 2 ) 

+ (C 3 + C A - C 6 )[Q b A(m 2 ,q 2 ) + Q s A(m 2 ,q 2 )]} . (3.5) 

Note that in this contribution the dominant charm-loop part has an enhanced Wilson co- 
efficient with respect to the one in the LO expression (3d). At q 2 < the hadronic matrix 



element A(m 2 ,q 2 ) is obtained from LCSR with .B-meson DA's, using the analytic expres- 
sion (4.8) in H derived for the c-quark case. The quark-flavour dependence is concentrated 
in the denominator which stems from the coefficient function of the effective operator. 

Finally, there is also a soft-gluon emission generated by the Og g operator as shown in 
Fig. ^|b. It is a nonfactorizable effect, so far not taken into account. In the next section 
we will derive a dedicated LCSR for the corresponding hadronic matrix element. 

Note that within the LCSR approach a systematic separation of "hard"- and "soft"- 
gluon contributions, with high and low average virtualities, respectively, is possible, at- 
tributing them to different terms in the OPE. E.g., in the sum rules with .B-meson DA's the 
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soft-gluon contributions enter the terms with three-particle DA's, whereas the hard-gluon 
contributions enter the NLO coefficient functions in the terms containing quark-antiquark 
DA's. Here we use LCSR's to compute the soft-gluon contributions whereas the hard-gluon 
effects are approximated by QCD factorization. 

3.4 Nonfactorizable spectator contributions 

A highly- virtual ("hard") gluon emitted from the intermediate quark loop or from the 
Os g operator vertex, can be absorbed by the spectator quark in the B — > K transition, 
generating an additional contribution - see Fig. |3|(b) or Fig. |||(c), respectively - which 
cannot be reduced to the B — )• K form factors. In || a factorization formula for this 
contribution was derived in terms of the hard-scattering kernels convoluted with the B- 
meson and kaon DA's: 

nj{BK) ( 2 n a S C F fBfKm b \ - f dlO ± f 1 (1) 

Konf'spect.il ) = 32nNc m | Z> J —^bH J q du {u)T ± (u, u) , (3.6) 

where the standard definition (see e.g., |l8|]) of two-particle B meson DA's <P^(oj) and of 
the twist-2 DA of the kaon ipx{u) is used. The NLO hard kernels T±\u,uj) corresponding 
to the sum of the diagrams in Fig. ||(b) and Fig. ||(c) are taken from || : 

T£\u,u>) = -— [Q c t||(n,m c )(C 1 +C4-C 6 ) + Q 6 t||(u,m fc )(C 3 + C4 
nib 

+Q s t ll (u,m s )(C 3 + C A - C fl ) + Q u t\\{uM.Ci - C 6 ) 
+Q d t ll (u,0)(C i -C 6 )] , 

T1'(u,uj) = -Q q k\ \g(m2,um% +uq z )(C 1 + C 4 + C 6 ) 

mew — q L y 6m\j L 

+g{m 2 b ,um 2 B + uq 2 )(C 3 + C A + C 6 ) 
+g{0, um% + uq 2 )(C 3 + 3C 4 + 3C 6 ) 

-^(C 3 - C 5 - 15C 6 )] + _^ 8 2/ 2 | , (3.8) 
3 u + uq A /m B J 

where is the electric charge of spectator quark in the B meson (q = u,d). The function 
tu(u,mq) convoluted with the S-meson DA 4>~^ arises from the two diagrams in Fig. ^(b) 
where the virtual photon is emitted from the quark loop. We recalculated and confirmed 
the resulting expression for t\\(u,m q ) presented in Eq. (28) of ||. In the remaining non- 
factorizable diagrams where the virtual photon is emitted from the quarks of the initial 
or final meson, only the DA <j) B contributes. Furthermore, we compared the above ex- 
pression with the dedicated sum rule calculation. To this end, a correlation function with 
.B-meson DA's and kaon interpolating current was introduced where the internal quark 
loop with the virtual photon emission was inserted. The resulting LCSR reproduces the 
QCD factorization expression ( |3.6D if two conditions are satisfied: (1) the asymptotic kaon 



-C 6 ) 

(3.7) 
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DA lpk(u) = 6u(l — u) is substituted in (|3,6|) and (2) in LCSR the leading-order two-point 
QCD sum rule for the kaon decay constant is used 2 . 

The factorization formula ( |3.6D by construction || only includes the leading power in 
the heavy-quark limit. We expect that LCSR can be more advantageous in assessing the 
power-suppressed corrections to this contribution. Note that in the factorization approach 
these corrections involve convolutions of the 5-meson DA's with higher-twist kaon DAs 
and generally suffer from end-point divergences. In the LCSR framework, albeit technically 
more involved, the off-shell correlation function is used where only the S-meson DAs enter. 

3.5 Weak annihilation 

For this contribution shown in Fig. |5] we use the same approximation as in || , employing 
the LO factorization formula 

<^tf> = 8j^pZ; / .hwMif («,«)■ (3.9) 

The hard-scattering kernels (u, oj) are given by 

rf (u,u) = o, 

Tl'(u,u) = e q ~ (C3 + 3C4). (3.10) 

mgu — g z nib 

The leading-power contribution - as well known - originates from the virtual photon 
emitted from the light-spectator quark inside the £>-meson. The other three diagrams in 
Fig. ||(a) are essential to ensure the e.m. current conservation. As noticed in [||], the 
weak annihilation diagrams develop an end-point singularity when the invariant mass of 
lepton pair becomes soft, q 2 ~ Aq CD . Hence working at large spacelike q 2 we avoid this 
problem. Radiative corrections to the weak annihilation diagrams (one of them shown 
in Fig. ||(b)) involving hard gluons are neglected since they have an additional 0(a s ) 
suppression. Note that a soft-gluon emission accompanied by weak annihilation apparently 
cannot be described by the above factorization formula but instead can be studied using 
LCSR's. Since there is yet additional power suppression involved we will not dwell on 
this problem here, expecting that altogether the above factorization formula provides a 
reasonable estimate of this suppressed effect. 

3.6 Power counting 

Concluding this section, let us discuss the power counting of separate contributions to the 
nonlocal amplitude T~L^ BK \q 2 ) considered above. As a reference scaling behaviour we take 
the nib ~^ 00 limit of the dominant factorizable contributions to B — > K£ + £~ amplitude 
Q2.2j ) which are determined by the B — >■ K form factors. Importantly, the large spacelike 
scale \q 2 \ 3> A.q CD introduced in our calculation does not scale with m&, hence it does 



2 Similar correspondence was found [[L9| between the LCSR for the B — > tt form factor in soft-collinear 
effective theory and the factorization formula for the hard-collinear spectator contribution to the same form 
factor. 
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not influence the familiar scaling behavior of the form factors in the large hadronic recoil 
region: fg K ~ fg K ~ 1/m^ 2 . This behaviour follows from the power counting in LCSR's 
with i?-meson DA's for these form factors, as well as from the alternative LCSRs with 
light meson (pion and kaon) DA's (see e.g.,[3C] for a detailed discussion). Note that in the 
numerical calculation here we use the B — > K form factors || from the LCSR's with kaon 
DA's calculated at finite m^, hence beyond the heavy-mass limit. 

Turning to T-L^ BK \q 2 ), we notice that the factorizable contributions (|3.1|) and (|3~g|) 



reveal the same power counting ~ 1/m^ 2 , since the coefficient functions stemming from 
the loops can only introduce a logarithmic dependence on the heavy mass scale. The scaling 
of the nonfactorizable soft-gluon contribution (3J3) is less trivial; it can be obtained from the 
heavy-quark limit of the LCSR for the hadronic matrix element A(m 2 , q 2 ) given in Eq.(4.8) 
in []|] . The derivation of the heavy-quark limit follows the same procedure as for the form 
factor LCSR's in Q. We obtain A(m 2 , q 2 ) ~ Aq CD /(Ai71 2 — q 2 ) x 0(m^ 3 ^ 2 ). Furthermore, 
the power-counting of the nonfactorizable spectator and weak annihilation contributions is 
directly obtained from the factorization formulae (3.6) and ( |3.9| ), respectively. Taking into 
account the scaling behaviour fg ~ l/m l J 2 and the fact that the kernels g and T^ 1,0 ^ in 
the leading power are independent of nib (modulo logarithms) we obtain that both these 
contributions scale as 0(m b ). Note however that numerically they are quite different, 



the weak annihilation contribution (|3.9| ) being suppressed due to the Wilson coefficients. 
Finally, the power-counting analysis of the soft-gluon nonfactorizable contribution of Os g 
presented in the next section reveals an additional 1/m 2 suppression with respect to the 
universal 1/m^ 2 scaling behaviour. Summarizing, the nonlocal amplitude 1-L^ BK \q 2 ) has 
the same leading-power behaviour as the dominant contributions of the 07^10 operators. 
The presence of the additional spacelike scale q 2 is crucial because the main nonperturbative 
correction to this amplitude stemming from the soft-gluon emission from the quark loops, is 
"protected" by ~ l/|g 2 | • The same Q 2 -scale plays a role of a regulator in the nonfactorizable 
spectator and weak annihilation contributions. This important circumstance justifies our 
approach of calculating the nonlocal amplitude at large spacelike q 2 . Future efforts to 
improve the adopted approximation, e.g., in the framework of LCSR's, should be invested 
in calculating the next-to-leading power corrections in 1/q 2 , such as the two soft-gluon 
emission diagrams. 



4. Soft-gluon contribution of the 0% g operator 

Here we present a calculation of the nonfactorizable effect in B — > K£ + £~ due to the 
soft-gluon emission generated by the gluon-penguin operator Os- The diagram in Fig.||(b) 
schematically represents this effect; important is that the gluon emission is a long-distance 
effect as opposed to diagrams Fig.^(a,c) where the gluon is absorbed at short distances and 
can be described by a perturbative propagator. The nonlocal hadronic matrix element 

i J d 4 xe^(K(p)\T{jUx),C e 8 ff O 89 (0)}\B(p + q)) = [(p • q)<f - 9 V]^gf (q 2 ) (4.1) 

includes all effects depicted in Fig. ||. As explained in the previous section, the soft-gluon 
part of this amplitude corresponds to the term in LCSR with the quark-antiquark-gluon 
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(c) 

Figure 7: Correlation function used to derive LCSR for the gluon-penguin operator contribution. 
The diagrams (a),(b) and (c) correspond to the hard-gluon and soft-gluon contributions, respectively. 



components of the B meson DA's. The derivation of the sum rule follows the procedure 
used in ||, and the operator Os g is simpler than the effective nonlocal operator of the 
soft-gluon emission from the quark loop. 

We start from the following correlation function 

U^(p,q)=i J d 4 y J d 4 xe^y + ^\0\T{j^(y),jU^,C e 8 ff Os 9 (0)}\B(p + q)) 

= ip a q^U(p 2 ,q 2 ) + ... , (4.2) 

where the gluon-penguin operator enters together with the e.m. current and kaon inter- 
polating current = d^y a ^^s. It is sufficient to consider the single kinematical structure 
shown on r.h.s., the rest is denoted by ellipses. At the hadronic level, we can express 
the above correlation function with the help of dispersion relation in the variable p 2 (the 
kaon-current momentum squared) at fixed q 2 : 

^(p, q) = l l KP ° \ {ip ■ qW ~ gV] «gf V) + ... , (4.3) 
m K ~ P 

where the ground-state term of the kaon contains the B —> K transition amplitude we are 
interested in, and ellipses denote the contributions from the excited and continuum states 
with the kaon quantum numbers. To derive the sum rule, we match the coefficient at the 
kinematical structure ~ p a q fl in the above to the invariant amplitude II(p 2 , q 2 ) in ( [4.21) . At 
sufficiently large \p 2 \ this correlation function ( f4.2| ) can be computed in terms of OPE in 
terms of .B-meson DA's. The gluon field emitted from the operator 0% g is either absorbed 
by one of the virtual quark lines, or enters the .B-meson state forming three-particle DA's 
(see Fig.|7|). Only the latter term in OPE interests us here, hence we only have to calculate 
the diagrams in Fig. ^(c). The diagrams in Figs. ^(a,b) correspond to the hard-gluon 
contributions already taken into account in terms of QCD factorization and included in 
the contributions 'r^f^^LO an< ^ %rwnfspect> as described in the previous section. 

We decompose the three-particle Fock state of the B meson in terms of four DA's: 
a(w , £) , ^V(u;,£), Xa(uj,£) and Ya(uj,£,), where u and £ are the momenta carried by the 
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light-quark and gluon inside the S-meson. The relevant definitions and ansatz for DA's 
(see App. B) are the same as the ones used in ||. After inserting this decomposition in 
the correlation function we calculate the sum of two diagrams in Fig. fj](c) with the virtual 
photon emission from b and s quarks; note that the photon emission from the spectator 
light quark does not contribute to LCSR. 

The invariant amplitude H(p 2 ,q 2 ) is obtained in the following form: 



[{p — uv)" 2 



n=l,2 ' 

F^ A ) (q 2 , U )9 A ( U , + FWO (q^^vico, £) 
-F^ (q 2 ,oj)X A (u } , + F^ (q 2 ,co)Y A (u;, £) 



(4.4) 



where Xa(w,£) = f (1i]Xa{i],Q, Ya.((jJ,£) = f df/Y^^, £). The explicit expressions of the 
o o 

coefficient functions Fn DA \q 2 ,uj) in = 1, 2) multiplying the three-particle S-meson DA's 

are collected in App. C. To simplify the calculation, we neglect the numerically very small 

~ £/m# terms in these functions, hence in the leading power they are independent of the 

gluon- momentum variable £. 

Equating two different representations of the correlation function, we obtain 

fK(m%-m 2 K -q 2 ) {BK) [°° p\s, q 2 ) , , 

In the above, the index "soft" at the amplitude 7i^ BK ^ indicates that we only take into 
account the soft-gluon contribution in the OPE result for the correlation function. 

The next important step is to use the quark-hadron duality approximation for the 
integral over the excited and continuum hadronic states with the kaon quantum numbers: 

t M.ir^|i, ,4.6) 

h S — V Z 7T Lk S — V 

where the effective threshold Sq is introduced. Finally, the Borel transformation is per- 
formed to improve the duality approximation and the final sum rule reads: 

rtfoSoM 2 ) = f , T yM J J ds ImU(s,q 2 ) . (4.7) 

J ' 8 fK{m z B -m z K -q 2 )ir J 

To obtain the explicit form of the integral on the r.h.s. it is convenient to start from the 



initial expression (4.4) and use the substitution rules presented in App. C for separate 
integrals in this expression. These substitutions simultaneously perform the transition to 
the dispersion form, subtraction of the higher states and Borel transformation. 

The soft-gluon contribution of Os g is suppressed by an extra power of 0(A 2 /m 2 ) with 
respect to the factorizable contributions, such as the LO loop contribution. Here A is 
a typical low-mass scale entering I?-meson DA's in HQET. The relevant power counting 
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can be carried out by expanding the sum rule (ft.7| ) at mj, — > oo. Focussing, e.g., on the 
contribution of the DA ^a(uj,^) in this sum rule, we obtain after using the substitution 
rules in App. C: 

.{BK)^ A ) ( 2 n m bfB r° du -e/M* 



"soft* * W 1 (ml - m 2 K - g2) ./„ 1 - u;/m B 



X 



/•oo 

/ d^F^ A \q 2 ,u) * A (u,,0 ~ , (4.8) 

•/ o 



taking into account that both the threshold parameter uq ~ Sq /mg and the coefficient 
function F-p^^ 2 ,^) scale as l/m&. 

5. Nonlocal hadronic effects in timelike region 

Collecting the results for separate contributions presented in the previous sections we cal- 
culate the nonlocal hadronic amplitude defined in ( |2.3| ) at q 2 < 0: 



To access the physical region Am 2 < q 2 < (ms — mjf) 2 of the B — > K£ + £~ decay, following 
§, we use hadronic dispersion relation in the channel of virtual photon. Here the following 
observation is important. The amplitude T~L^ BK \q 2 ) in (5.1) can be decomposed in separate 



contributions distinguished by the flavour of the quark interacting with the virtual photon 
via e.m. current. One simply collects all terms in H^ BK \q 2 ) that are proportional to 
the electric charge Q q of the quark q. This flavour splitting allows one to establish a 
correspondence between the part of ( |5.1| ) with a given Q q and the hadronic states with 
the quark content qq in the dispersion relation. For example, the part of the amplitude 
H( BK \q 2 ) proportional to Q u ,d (Qs) is dual to the part of the hadronic dispersion relation 
containing the vector mesons p and oj ((f)) and excited and continuum states with the 
same quantum numbers. Here we neglect small mixing effects. Accordingly, the part of 
the amplitude 1-L^ BK \q 2 ) proportional to Q c generates the vector charmonium states and 
open charm states in the timelike region. Since q 2 < m\ in B — > K£ + £~, the part of 
the amplitude proportional to Qb, only contributes to the nonresonance background in the 
dispersion relation. 

We separate the contributions from different flavors of quarks in ( |5.1[ ) in the following 

way: 

%^tf)=uf«\ q 2 ) + U% K \ q 2 ), (5.2) 

so that the first (second) term in the above contains all contributions from the different 
terms in ( |5.1[ ) proportional to Q s , c ,b (Qu,d)- The reason why the s,c, 6-quark contribu- 
tions are not separated from each other is that in certain terms included in ( |5.l[) e.g., in 



^Mrf atlo(^ 2 )' * ne vhtual photon emissions from b and s quarks are interconnected by 



gauge invariance. Also the results of [17] used to calculate these contributions do not allow 



one the separation of the c-quark component. In future, a more detailed flavour separation 
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can be achieved, performing a dedicated analysis of the NLO two-loop diagrams. Another 
important point is that the photon emission from the spectator light quarks also contributes 
to T~L^ BK ^ (q 2 ). Hence, one has to specify the light flavour of the decaying B meson. Our 
default choice is B° -> K°£+£-. To calculate the isospin asymmetry we will also consider 
B~ —7- K~l + l~ . In the approximation adopted in this paper, the small CP violation effects 
are neglected, hence the same results are valid for B° and B + decays, respectively. 
The separate hadronic dispersion relations for the two parts of the decomposition 

are: 



^ BK )( f fi\-'lji BK )(n 2 \^(n 2 a 2 \\ V Kyf V \A BVK \e l 
H scb (9b) + (9 -9b) 2^ (m 2 2 )(m2 2 _ • ™ 



(s - gg)(s - <? 2 - ie) 



(5.3) 



and 

,(£*0/ 2\_nj(BK) ( 2 n . / 2 2n T /y | ABVid eW 



, (9 2 ) = ^ ) (9o 2 ) + (9 2 -9o 2 )[ £ 



P, \r ds - 



S (s-9o)('S-9 2 -«e) 



(5.4) 



In the above, the decay constant of a vector meson 1/ with the polarization vector ey is 
defined as: (0\jem\V) = Kymy fye v \ the coefficients k p = l/v2, k w = l/(3\/2), = —1/3, 
K J/ip = K i>{2S) = 2/3 follow from the valence quark content of V. As discussed in Q 
one subtraction at q$ guarantees the convergence of the dispersion integrals. We choose 
<7q = —1.0 GeV 2 and use the calculated result to fix the subtraction terms 'H^buAlo)- 
The advantage of the flavour separation done above is that the suppressed part of the 
nonlocal amplitude, H^4 K \q 2 < 0), has its own dispersion relation, allowing one a more 
accurate estimate of contributions with p, oj quantum numbers in the timelike region, which 
otherwise would have been hidden under the dominant contributions of charmonium states. 

The residues of vector meson poles in both dispersion relations ( |5.3D and fl5.4j ) contain 
amplitudes Abvk of the nonleptonic decays B — > VK. Final-state strong rescattering is 
taken into account by the phases attributed to each amplitude 3 . The B — > VK decays 
are well measured, allowing one to extract the absolute values of their amplitudes. The 
relative phases of these amplitudes are unknown and they may considerably influence the 
pattern of vector meson contributions in the dispersion relations. Note that the absolute 
values Ivlwid extracted from the experiment deviate from their naive factorization esti- 
mates obtained multiplying fy with the combination of Wilson coefficients and the B — >• K 
form factor. These deviations are considerable for V = J/ip,ip(2S), revealing large non- 
factorizable effects in B — > J /tpK,tp(2S)K , they are somewhat smaller for the decays into 



3 As discussed above, these phases are related to the discontinuities in the second variable (p + q) 2 . 
In this sense the dispersion relations we are considering here can be interpreted as a double-dispersion 
representations for the amplitude H' SK ' ((p + q) 2 , q 2 ), where each amplitude Abvk is itself expressable in 
a form of single dispersion relation over (p + q) 2 . 
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light vector mesons. Altogether, the dispersion relations (|5.3| ) and (|5.4j ) contain too many 
hadronic degrees of freedom to be fixed and/or constrained solely from the matching of 
the r.h.s. with the result of the calculation at negative q 2 . Therefore, following the ap- 
proach advocated in Q we fix the absolute values of residues in both dispersion relations, 
employing experimental data on the decay constants and nonleptonic decay amplitudes. 
The main difference between the relations considered in P] and here is that, in addition to 
the contributions of charmonium states, we also include the contributions of light vector 
mesons employing the data onB-> (pK and B -4 p(u)K partial widths. Importantly, the 



data currently accumulated in |20|, |21[ allow one to fix the decay amplitudes for both B" 
and B~ mesons separately, which is particularly important for the evaluation of the isospin 
asymmetry. 

After fixing the residues of the lowest resonance poles in both dispersion relations (|5.3| ) 
and ( |5.4[ ), they still contain a significant amount of unknown hadronic degrees of freedom 
from higher states, accumulated in the integrals over spectral densities. Note that the 



integral in (5.3) contains the contributions of excited c/>-resonances and continuum states 
with the same quantum numbers, as well as the contributions of all vector charmonium 
states and open charm-anticharm states above the threshold s = Am 2 D . The integral in ( p.-j| ) 
in its turn, accumulates excited and continuum states with the p, uj quantum numbers. At 
very large s, approaching the B* meson mass squared the bs states contribute in both 
integrals. At low and intermediate q 2 the presence of these states is signaled (at least in 
the factorizable part of Ti^^^q 2 )) by the growth of the B — > K form factor. The fit of 
the dispersion relation to the calculation result at negative q 2 will introduce important 
constraints on the hadronic integrals entering dispersion relations. However, in order to 
perform a continuation to positive q 2 one needs to parametrize these integrals in terms of 
a certain ansatz. 

Let us first concentrate on the dominant part of the amplitude parametrized by the 
hadronic dispersion relation (j5.3|). Note that the region of q 2 we are interested in, is located 
below the open charm threshold Amjj. The only (subdominant) contribution in this region 
to the spectral density p sc b{s) originates from the superposition of ss states starting from 
(f> meson and including its excitations and multiparticle states, such as KK with J p = 1~. 



The ansatz used in QCD sum rules [22] to describe the ss -channel with vector quantum 
numbers is the <j) -meson pole combined with the continuum, the latter approximated by 
the quark-parton duality with an effective threshold. This motivates us to replace the 
spectral density by the following ansatz: 

Psc b (s)6(s - 4) = -hrmfi(s)e(s - 8+) , (5.5) 

7T 

where the effective threshold Sq is determined from the QCD sum rule in the cj) channel. 
The simplest choice for the r.h.s. is to approximate it with the LO factorizable part 
U-fact lo s wri ere the index s indicates that only the part proportional to Q s is taken in 
( |3.1| ). The spectral density is then reduced to the imaginary part of the s-quark loop 



function (^2) multiplied by the B — )■ K form factor. Note that in this case we adopt a 



stronger assumption of local duality, rather than the semi-local duality used in QCD sum 
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rules. Also at s < Am 2 D , that is within the validity region for LCSR, the calculated form 
factor can still be used. We therefore subdivide the integral in the dispersion relation (|5.3| ) 
in two parts and adopt the approximation for the interval Sq < s < Am 2 D . In addition, 
the imaginary part of the integral at q 2 > Sq is modified by introducing an effective width 
factor in the denominator, replacing s — q 2 — it — >■ s — q 2 — iy/s T e ff(s). The dependence of 
this width factor on s is taken as T e ff(s) = ^^fs®(s — Arn 2 K ), following the resonance model 
adopted for the timelike form factors in |23| where the universal value 7 = 0.2 was derived 
from the spectrum of light vector mesons. The step-function accounts for the kinematical 
limit where the width vanishes. 

In the remaining integral from 4m 2 D to infinity, the spectral density p sc b{s) cannot 
be simply parametrized because it contains a complicated superposition of four hadronic 
components: 1) broad charmonium resonances and open-charm continuum states (the 
dominant contribution), 2) the "tail" of ss states, 3) the states with B* quantum numbers 
and 4) the bb states corresponding to the photon emission from b quark. Note however, that 
in the region of our interest q 2 < Am 2 D (practically, below the J/ip mass), this integral has 
no singularities and can be represented as a generic series expansion in powers of q 2 /Am 2 D . 
An alternative is to perform a usual q 2 — >■ z transformation and to use the z-expansion. 
As a default model we use a rather simple approximation 

f°° ,1 P scb ^ j_ h J? (kr\ 

/ as- 2T7 5 rr- = a S cb + b sch - — T , (5.6) 

where a sc b and b sc b are unknown complex constants. 

Turning to the second dispersion relation (|5.4|), we use a similar approximation for the 
spectral density: 

Pud(s)9(s - S§) = hvM*?{s)e{s - ag) , (5.7) 

where Sq is the effective threshold, typical for the QCD sum rules in the p channel. Ac- 
cordingly, the integral over large s is represented similar to (|5.6| ) with the two additional 
complex parameters a u d and b u d- For the sake of uniformity we put the lower limit of this 
integral to the same value Am 2 -). 

After parametrizing the integrals over the spectral densities, the dispersion relations 
fl5.3| ) and ( |5.4| ) are fitted to the results for T-Lf c f(q 2 ) and H^d (^ 2 ) respectively, at negative 
q 2 . The fitted parameters are the coefficients a<j C ;,,5 sc b in the expansion of the integral 
Q5.6| ) and a U( [, b u d for the corresponding integral for u, d quark states. After that, the 
two dispersion relations are added together and continued to timelike q 2 resulting in the 
desired nonlocal amplitude % BK {q 2 ) for B — > K£ + £~ decay. To assess the dependence 
on the parametrization of the dispersion relations, apart from the default parametrization 
of the dispersion integrals described above, we also considered three alternative versions: 
(I) the integral ( |5.6|) is replaced with an effective pole as it was done in ||; (II) a generic 
z-series parametrization for the integral over higher states (|5.6| ) is used, and (III) as an 
extreme choice only the dominant J/ip and ip(2S) states are left in the dispersion relation 
with no flavour-splitting applied. 
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Figure 8: Nonlocal hadronic amplitude in B — > K£ + £~ , expressed in terms of the correction 
to Cg, calculated at q 2 < {or the central values of the input. The separate contributions are 
from the factorizable quark loops (ACgf, black), weak annihilation (ACgwA, dashed black), soft- 
gluon nonfactorizable emission (ACg n f, green), soft-gluon emission from Og g operator (ACgc, 
multiplied by 10 2 , orange), real part of NLO hard-gluon effects (Re[ACg nlo], brown), total real 
part (Re[ACgtot\, red), imaginary part of NLO hard-gluon effects (Im[ACg nlo], violet). 



6. Numerical results 

We start with the numerical evaluation of the amplitude ( |5.1[ ) in the region q 2 < 0. The 
same input as in Q is adopted, where one can find more details. In particular, since the 
calculation is done at spacelike momentum transfer, it is more appropriate to use the MS 
quark masses, for which we adopt the following intervals: m s {2 GeV) = (98 ± 16) MeV, 
m c (m c ) = (1.29 ± 0.03) GeV and m b {m b ) = (4.164 ± 0.025) GeV. The functional form of 
-B-meson DA's entering the factorization formulae and LCSR's is specified in App. B. We 
use the interval f B = 180 ± 30 MeV for the decay constant, A#(l GeV) = 460 ± 110 MeV 
for the inverse moment and Ag(l GeV) = y/3/2 Ab(1 GeV) for the parameter entering the 
3-particle DA's. For the kaon decay constant we use Q f K = 159.8 ± 1.4 ± 0.44 MeV 
and for the Gegenbauer moments of the kaon DA: af"(l GeV) = 0.10 ±0.04, a^(l GeV) = 
0.25 ±0.15. Finally, in LCSR with the kaon interpolating current the Borel parameter and 
effective duality threshold are taken as M 2 = 1.0 ± 0.25 GeV 2 and s$ = 1.05 GeV 2 , 
respectively. The default value of the renormalization scale is fj, = 3.0 GeV, so that 
a s (3 GeV) = 0.252, and the scale is then varied within 2.5 GeV < fx < 4.5 GeV. 

For the B — )■ K form factors the predictions from LCSR with kaon DA's are used. 
They are presented in the App. B of ||, in particular, the normalization of the vector and 
tensor form factors are given by fg K (0) = 0.34^q'q 2 , and f]^ K {0) = 0.39^003, respectively. 
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In the contributions calculated using QCD factorization approach we adopt a universal 
B — ^ K form factor £,bk(q 2 ) which is set equal to /^(<? 2 ), neglecting the small difference 
between tensor and vector form factors. 

According to (|2.4|), the result for T~L^ BK \q 2 < 0) is normalized to the form factor, 
yielding a correction to the Wilson coefficient Cg. Separate contributions to this correction 
are plotted in Fig. |8| for the central values of the input. Numerically, in the real part 
of AC { g BK \q 2 ) there are substantial cancellations between separate contributions, and 
the soft-gluon nonfactorizable part plays an important role as already observed in ||. 
Altogether, the real and imaginary parts of ACg BK \q 2 < 0) - the latter generated by the 
NLO hard-gluon effects - reach the level of a few percent of the short-distance coefficient 
C 9 . 

To assess the impact of nonlocal hadronic effects on the observables in B — > K£ + l~, 
we turn now to the numerical analysis of the dispersion relations, allowing us to smoothly 
continue these effects into the physical region q 2 > 0. To this end, we follow the procedure 
described in the previous section, splitting 7-L^ BK \q 2 ) according to (fT2]), and fitting the 
functions Ti^^iq 2 ) and H^^iq 2 ) in the region —4m 2 < q 2 < —1.0 GeV 2 to the two 
dispersion relations ( |5.3| ) and ( |5.4| ), respectively. In these relations, the absolute values of 
the nonleptonic B — > VK amplitudes, together with the decay constants of vector mesons 
entering the residues of the resonance poles, are extracted from the experimental data [20] 
on the B — > VK and V — > widths. These parameters are collected in Table |]. 
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Table 1: Decay constants of vector mesons and amplitudes of B 
calculated from the experimental data \2L 



VK decays (all in MeV) 



Furthermore, the integrals over the spectral functions of higher states in (5.3) and (|5.4| ) 
are subdivided in two parts. The integral below Am 2 D is parametrized employing the local 
duality approximation with the effective thresholds Sq = 1.95 GeV 2 in ( |5.5| ) and Sq = 1.5 
GeV 2 in (5.7), whereas for the remaining integrals above 4m 2 D the polynomial approxima- 
tions in the form ( |5.6| ) are used. In ( |5.3| ) the phases of J/ij) and i^(2S) contributions are 
varied independently at — tt < ^Pj/^,f^(2S) ^ n an d the 4> contribution is taken real. For 
each combination of the phases (pj/-^ and ipM^s) the fit of the complex parameters a sc b 
and b sc b is repeated. The best fit is obtained when a destructive interference between the 
J/ij) and '0(25') terms in ( |5.3| ) takes place, as observed before in |^]. The second dispersion 
relation (|5.4j) is treated in a similar way. The phases and fitted parameters are collected in 
Table |2[ Substituting them in the dispersion relations, and continuing the latter to q 2 > 
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Figure 9: The effective correction ACg BK \q 2 ) in the physical region of B — > K£ + £~ : the red (blue) 
solid curve corresponds to the real (imaginary) part obtained from the hadronic dispersion relation, 
fitted to the QCD calculation at q 2 < (central input, default parametrization) . The shaded areas 
indicate the uncertainties. The dashed curves correspond to the prediction of QCD factorization 
obtained with same input. 



we finally obtain a numerical result for the nonlocal hadronic amplitude rl BK {q 2 ) in the 
physical region. 

The resulting effective correction AC BK (q 2 > 0) is plotted in Fig. || in the region of 



large hadronic recoil, up to q 



~ m 



At q 2 < m 2 , the behavior of AC BK (q 2 ) reflects the 



presence of light vector resonances. In the same figure we compare our result with the result 
for AC BK directly calculated at q 2 > 0, using the same approach and approximation as 
in p]. This comparison is possible in a restricted region, 2.0 GeV 2 < q 2 < 4.0 GeV 2 , above 
the light resonances and sufficiently below the cc-quark threshold. In this region, the real 
parts AC B K (q 2 ) obtained from our dispersion approach and from the direct calculation are 
in agreement within uncertainties, whereas the imaginary parts deviate from each other. 

The uncertainties shown in Fig. |9] are estimated by varying each input parameter 
involved in the calculation of /S.Cq K (q 2 < 0) and then recalculating ACg K (q 2 > 0) from 
the dispersion relation with the modified input. The resulting individual deviations of 
ACg K (q 2 > 0) are then added in quadrature. Furthermore, we repeat the fit of the 
dispersion relations for the three alternative parametrizations of the dispersion integrals 







\a scb \ [GeV" 2 ] 


Arg[a scb ] 


|6 sc6 |[GeV- 2 ] 


Arg[6 scf) ] 


-2.14 


0.77 


1.40 x 10" 5 


0.69 


1.37 x 10" 4 


0.86 






M [GeV" 2 ] 


Arg[aw] 


Kd\ [GeV" 2 ] 


Arg[6 ud ] 


0.64 


-2.50 


2.74 x 10~ 5 


-2.13 


5.75 x 10- 5 


-2.07 



Table 2: Parameters of the dispersion relations (5.3) and (5.4) obtained from the fit to the ampli- 
tudes 'rl[cb < ^ an( l ^ud^ at g 2 < for the central values of the input. 
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Table 3: The effective correction to the coefficient Cg at different q 2 . 



n 2 (CpV 2 \ 


Re.r /\r BK t n 2 w 

_ie[_vO g [q )\ 




1.0 


n 1 1 +°- 06 


n 33+0.11 


9 n 
z.u 


D 1 rs+0.07 
u - lo -0.06 




9 c< 
_.o 


17+0-07 
U - L '-0.07 


n qn+o-ii 
-u._j_ 013 


O.U 


fl 1 8+0-08 
-U.18_ 008 


n Q1+0-12 
-U.6L_ 0U 


O.U 


n 9f) +0.09 


n Q9+0-12 
u - oz -0.16 




n 91 +0.09 

-u._i_ 12 


fl 34+0-12 


-.0 


n on+0.10 
u - zz -0.13 


fl 3R+0-13 


O.U 


n 93+0- 11 

.15 


n 3Q+0-13 

u.oy_ 22 


o.o 


o ?^+°- 12 

u - zo -0.17 


n 42+0-1 3 

u - 4z -0.25 


6.0 


-0.281°;^ 


-0.48^ 


6.5 


-0.32^- 


-0 ^+°-l 5 
u - oo -0.30 


7.0 


-0-38±-it 


-0.66±g5i 


7.5 


-0.48™* 


-0.83_™ 


8.0 


-0 65+ - 18 


1 19+0.22 

1 - lz -0.40 


8.5 


-1 0f +0 - 21 


-1 6Q+ - 31 
j-.oy_0.45 


9.0 


-2 02+ - 34 
z - uz -0.38 


_3 2Q+0- 57 
°- zy -0.56 



(I)- (III) specified at the end of the previous section. We interpret the difference between 
ACg K (q 2 > 0) obtained with the default parametrization and the one in the model III 
(which only contains J /tp and ip{2S) resonances) as the "systematic" uncertainty of our 
approach and include it in the uncertainty budget. As an example, we present the correction 
to Cg calculated at one particular value of dilepton mass: 



Re[AC^ K {q 2 = 4GeV 2 ) 



ImlACfV 



4 GeV" 



-0.21; 



-0.34 



+0.03 




fO.OO 




f0.08 




f0.02 




f0.02 


-0.04 


As 


-0.10 


af 


-0.00 


I-' 


-0.01 


/ + x(0) 


-0.06 



syst. 



(6.1) 



-0.05 
-0.07 



+0.01 




+0.08 




+0.06 


+0.05 




+0.02 


-0.12 




-0.07 




-0.03 


-0.02 




-0.04 


As 


af 




af 






/sk(0) 





syst. 



where all significant (> ±5%) uncertainties are shown separately. As expected, the de- 
pendence on the B — )• K form factor is inessential. Here, we assume a fixed uncertainty 
for B — > K form factors stemming from LCSR and neglect possible correlations when 
varying the rest of the input. For convenience, we also present the numerical intervals for 
ACg K (q 2 ) in the region 1.5 GeV 2 < q 2 < 9.0 GeV 2 in Table |. They have to be compared 
with Cg ~ 4.0 — 4.5 (see App. A). The magnitude of ACg K (q 2 ) remains in the ballpark 
of ~ 10% of Cg at small and intermediate q 2 and grows approaching the charmonium re- 
gion. We refrain from quoting ACg K (q 2 ) in the region between J ftp and ipi^S) where the 
"systematic" uncertainty caused by the parameterization of dispersion relation is rather 
large. 
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7. Observables in B -+ K£+t 



we calculate the 



K°n + fi in the region 4m^ < q 2 < Am 2 D (see Fig. [lQj ) . Integrat- 



Substituting the numerical results for T-L^ BK \q 2 ) to the amplitude 
differential width of B° 
ing over the typical intervals of dilepton-mass squared (bins) selected in the experiments, 
we present the resulting partial widths in Table |4], in comparison with the available mea- 
surements. In order to illustrate the origin of the quoted uncertainties, we present them 
separately for one of the bins: 



;.o GeV 2 



dq 



2 dBR(B° -> K°n + fi- 



dq 2 



1.0 GeV 



1.76 



+0.58 


+0.16 


+0.01 


+0.01 




-0.21 


-0.09 


-0.01 


-0.02 


syst. ^ 




fi K (°) 


slope 


H 





x 10 



-7 



(7.1) 



where the remaining individual uncertainties smaller than 1% are not shown. As opposed to 
ACg BK \ the theoretical uncertainty of the width mainly originates from the normalization 
of the form factor fg K - Our predictions for the partial widths of B° -+ K°£ + £~ in Table [| 
are somewhat larger than the results of the most recent LHCb measurement |Q] of B~ — > 
K~£ + £~. 



Table 4: dBR(B — > K^ + fi )/dq 2 integrated over [q^iri-i Qmax] lTl units of 10 7 . 



\flmin i Qvnax\ 

GeV 2 


Belle § 


CDF § 


LHCb § 


LHCb § 


this work 


[0.05,2.0] 


0.81±g;i| ± 0.05 


0.33 ±0.10 ±0.02 


21 + 027 

u - zi -0.23 


0.56 ±0.05 ±0.03 


o.7i±g:g 


[2.0,4.3] 


0.46±g;^ ± 0.03 


0.77 ±0.14 ±0.05 


n n7 +0 - 25 

u - u ' -0.21 


0.57 ±0.05 ±0.02 


0.80t°; 27 


[4.3,8.68] 


1.00±g;£| ± 0.06 


1.05 ±0.17 ±0.07 


1.2 ±0.3 


1.00 ±0.07 ±0.04 


1 QO+0.53 
i.O»_Q 2 2 


[1.0,6.0] 


1.36±g:lf ± 0.08 


1.29 ± 0.18 ± 0.08 


0.65±«;|1 


1.21 ±0.09 ±0.07 


1 7fi+ - 60 
i - '°-0.23 



In addition, we also calculated the (CP averaged) isospin asymmetry in B — > K£ + £ 
defined as 

(0-) , 2) = dr(B -» Ko£+£-)/dq 2 - dT{B- -> K-£+£-)/dq 2 
G/ (q ' dF(B -»• K £+£-)/dq 2 + dT{B~ -> K~£+£-)/dq 2 ' 1 ' ' 

In our approach this effect originates from the small differences between the amplitudes 
of photon emission from the d and u spectator quarks, in the nonfactorizable NLO con- 
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2 4 6 8 10 12 14 

q 2 (GeV 2 ) 

Figure 10: Differential branching fraction of B° — > K°[i + . The solid line corresponds to the 
central input and to the default parametrization for the dispersion integrals. The darker (green) 
and brighter (yellow) shaded area indicates the uncertainties including (excluding) the one from the 
form factor normalization. The alternative parametrizations I, II and III of the dispersion integrals 
yield dotted, dash-dotted and thin dash- double- dotted (blue) lines, respectively. The long-dashed line 
corresponds to the width calculated without nonlocal hadronic effects. 



Table 5: Isospin asymmetry calculated from the partial widths of B° — > K°£ + £ and B 
K~£+£- integrated over 1.0 < q 2 < 6.0 GeV 2 . 



Belle § 


BaBar || 


LHCb § 


this work 


— 0.41±g;|g ± 0.07 


-0.41 ±0.25 ±0.01 


-o ^+ - 23 


(-0.4)% + (-0.3)% 



tributions and in the weak annihilation. The result for the differential asymmetry defined 
above is shown in Fig. 11 and does not exceed a ±4% level within estimated errors. At 
q 2 > 2 GeV 2 our expectation for the isospin asymmetry in B — > K£ + £~ is in the ballpark 
of the prediction obtained in [24] for the isospin asymmetry in B — > K*£ + l~ . In Table ||| we 
present the integrated isospin asymmetry defined similarly to (|7.2j) but with the partially 



integrated widths instead of differential widths. This integrated characteristics is expected 
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Figure 11: Isospin asymmetry af \q 2 ) in B — > K£ + £ (the solid curve), the shaded (green) band 
indicates the uncertainty of our prediction. 



to be very small, revealing an intriguing tension with the available experimental measure- 
ments presented in the same table. It is important to further improve the accuracy of our 
prediction. E.g., the calculation presented here assumes isospin symmetry for the B — )■ K 
form factors. This issue deserves a separate study but - having in mind the usual mag- 
nitude of isospin violation - it is hard to expect that the deviation of fiojso from fi-v- 
brings substantial changes in the isospin asymmetry of B — > K£ + £~. 

Having at hand the SM prediction for the amplitude of B — > K£ + £~ it is interesting 
to investigate the potential influence of new physics in b — > s£ + £~ on the observables in 
this decay. Recent improvement [25] of the upper bound for BR(B S — > /i + ^~) has already 



put substantial constraints on the new FCNC operators involving pseudoscalar and scalar 
couplings of the lepton pair. Hence we only consider a generic new physics induced by the 
two tensor operators: 



a. 



2tt 



T5 



Or 



2tt 



(7.3) 



with Wilson coefficients Ct and Ct5, respectively. 

The B —7- K£ + £~ width is calculated adding the new operator contributions to the 
decay amplitude and assuming Ct = ±Ct5- As shown in |1C], these two parameters are 



effectively constrained by comparing the measured upper bounds on the inclusive branching 
fraction B — > X s £ + £~ with the SM prediction. Typically one obtains \Ct\ ~ |Qrs| < 1>2. 
With this constraint, our result for the decay rate of B — > K£ + £~, including the new 
physics contribution and integrated over the interval 1.0 < q 2 < 6.0 GeV 2 , reveals a rather 
small deviation from the SM prediction, at the level of < 5%. The predicted forward- 
backward asymmetry in B — > K£ + £~, emerging due to the new operators (7.3), is plotted 
in Fig. |l2] at Ct = Ct5 = 1.2. This observable can reach at most —(5 — 10)%, being 
almost independent of q 2 . There are measurements of the forward-backward asymmetry 
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Figure 12: Forward-backward asymmetry in B — > Kt + £~ emerging due to new physics contribu- 
tions of tensor operators with the Wilson coefficients Ct = C*T5 = 1.2. The shaded (green) band 
indicates the uncertainty of our prediction due to variation of SM parameters. The data points are 
from the measurement by LHCb B (triangles) and Belle (circles). 



in B — > K£ + £ by Belle Collaboration || and more recently by LHCb Collaboration Q, 
see Fig. [l^. However, these measurements are not yet sensitive to such small effects. 

8. Discussion 

In this paper, we calculated the nonlocal hadronic contributions to B — > K£ + £~ generated 
by the four-quark and gluon-penguin operators of the effective Hamiltonian, combined 
with the e.m. emission of the lepton pair. In this calculation the LCSR results for the 
form factors and for the soft-gluon nonfactorizable effects are used, whereas the hard- 
gluon effects are approximated employing QCD factorization. We followed the method 
suggested in combining the QCD results for the nonlocal hadronic matrix elements 
valid in the region of negative dilepton-mass squared, with the dispersion relations in the 
physical region. In these relations the residues of the lowest vector mesons are fixed from 
the data, whereas a nontrivial destructive interference between the dominant J/tp and 
^(25) contributions plays an important role. A model ansatz is adopted for the higher 
state contributions to the dispersion relations, with the parameters fitted to the nonlocal 
amplitude calculated in QCD. 

Our main result is the numerical prediction for the nonlocal hadronic amplitude 
cast in the form of the correction to the Wilson coefficient Cg. This correc- 
tion displayed in Table ^ and in Fig. |9| can be used in future phenomenological analysis 
of the B —7- K£ + £~ decay. Employing dispersion relation, we avoid unphysical "kinks" in 
the observables generated by the quark-antiquark production thresholds which appear if 
one directly uses QCD diagrams in the physical region (see e.g., [^6|). Probing alternative 
parametrizations of the dispersion integrals we argue that the "systematic" uncertainty 
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introduced by this procedure is inessential in the large hadronic recoil region, practically 



up to q 2 = m 2 j, , . The approach used here is not applicable at larger values of dilepton 



mass, especially at q 2 > irn/jj, where a complicated interference of charmonium resonances 
and the proximity of the singularities of the B —> K form factors brings too many hadronic 
degrees of freedom into play 4 . 

In the large hadronic recoil region, we evaluated the partial width and isospin asym- 
metry of B — > K£ + £~ . The impact of nonlocal hadronic contributions on the width is 
moderate below the charmonium region. The main uncertainty of our prediction in this 
region stems from the normalization of the B — > K form factors, leaving a room for im- 
proving the accuracy of the SM prediction for B -4 K£ + £~ . In particular, new precise 
lattice calculations of these form factors combined with updated LCSR results at small 
q 2 are desirable. In addition, we predict a very small isospin asymmetry in B — > K£ + £~ , 
hence, more precise measurements of this observable are needed. 

A further improvement of our prediction for the effective correction ACg BK \q 2 ) is 
possible along two lines: firstly, achieving more accuracy in QCD calculations at negative 
q 2 and secondly, with more detailed parametrizations of the dispersion integrals adding the 
excited vector meson terms. In this respect, dedicated measurements of nonleptonic decays 
of the type B —> V'K where V 1 are radial excitations of p,uj,<p and the charmonium levels 
starting from -0(3772) will be helpful. 

The study presented in this paper could be extended to B — >■ K*£ + £~, for which the 
charm-loop effect has already been analyzed in ||. However, calculating the other effects in 
this process will be more demanding. E.g., the nonlocal hadronic contributions proportional 
to Q s such as the ^>-meson pole will be enhanced at small q 2 , being multiplied by the factor 
1/q 2 from the virtual photon propagator in the amplitude for the transversely polarized 
K* meson. To accurately analyze this particular contribution, a separation of Q s — Qb 
component from the Q c -one is desirable in all three nonlocal hadronic amplitudes Ti\ 2 3 
corresponding to the three polarization states of the K* meson. This demands a dedicated 
"flavour splitting" of the NLO diagrams calculated in [16, [l?]] on one hand and a separate 



dispersion relation with ss vector-meson states on the other hand. Also the NLO two-loop 
diagrams induced by the penguin operators should be included to achieve an adequate 
accuracy. Altogether, the QCD calculation at negative q 2 and a complete analysis of 
hadronic dispersion relations for all three invariant amplitudes of B — > K*£ + £~ represents 
a significantly more challenging task. The nonlocal corrections to Cg are expected to be 
numerically larger than for the kaon mode, as already seen from the comparison of charm 
loop effects for both processes in ||. Hence, a substantial impact of the nonlocal hadronic 
effects on B — > K*£ + £~ is anticipated, even in the large hadronic recoil region. Last but 
not least, as already discussed in the introduction, the fact that K* has a finite width 
decaying to Kir, hinders one from the calculation of "pure" B — > K* form factors, with 



4 Let us mention that in the low hadronic-recoil region a different method is used |27j], based on the 
local OPE of the nonlocal hadronic amplitudes valid at q 2 ~ mj; — > oo. Applying this method to B — > 
K^k+r decay, one relies on the (quasi)local quark-hadron duality and smallness of power corrections. 
The connection of the local OPE to the approach used in this paper is an interesting open problem that 
deserves a dedicated study. 
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the accuracy achieved for the B — >■ K form factors, either on the lattice or with QCD sum 
rules. New approaches are desirable to calculate the form factors of S-meson transitions 
to the two-meson system including resonances. Only in this case the theory can meet 
the challenge of continuously improving experimental measurements of exclusive FCNC B 
decays. 
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Appendix A: Effective Hamiltonian 

Here we list the operators entering the effective Hamiltonian ([2.1|): 



0\ = (sli p cl) (cl7 p 6l) 
5 = (sLiph) E (qri P( ir) 



o 7 , 



Y^sa^inisL + m b R)bF' 
O9 = ^ (s Llp b L ) (Pfl) 



47T 



s l1p c l 



2 



do 



16tH 



a,. 



4-7T 



(sLl P b L ) {h p 1sl) , 



where the notations are qum = 1- ^ 75 q and L(R) = 1 ~^^ 75 . We use the standard 
conventions for the operators Oi, except the labeling of 0\ and O2 is interchanged. The 
sign convention for Oj^ and Os g corresponds to the covariant derivative iD^ = id^ + 
eQfA^ + gT a A^, where Qj is the fermion charge. In addition, the convention for the 
Levi-Civita tensor adopted in this work is Tr^YY^J 5 } = 4ie^ pA , e 0123 = — 1. The 
numerical values of the Wilson coefficients at three different values of the scale [i are listed 



in the table below [13]: 



H (GeV) 


0.5rrib 


m b 


1.5m;, 


Ci 


1.193 


1.117 


1.090 


c 2 


-0.401 


-0.267 


-0.214 


C7 3 (xl0- 2 ) 


1.919 


1.206 


0.931 


C 4 (xlO- 2 ) 


-3.964 


-2.750 


-2.225 


C 5 (xl0- 2 ) 


1.041 


0.770 


0.639 


C 6 (xl0" 2 ) 


-5.479 


-3.411 


-2.637 


c 7 


-0.370 


-0.320 


-0.297 


C 8 g 


-0.184 


-0.166 


-0.157 


c 9 


4.450 


4.232 


4.029 


C10 


-4.410 


-4.410 


-4.410 
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We also use \V ts \ = 0.04051™ , \V tb \ 



O-QSQiaooooi HI and a em = 1/129. 



Appendix B: B -Meson Light-cone DA's 

In x 2 = limit we adopt the following decomposition of the l?-to-vacuum matrix element 
into four independent three-particle DA's (see e.g., [29(]): 



(O\q a (x)G Xp (ux)h v0 (O)\B°(v)} 



fBm B 



oo oo 



duj J d£ e 
11 o 



-i(uj+u^)v-x 



\ f • x 



V ■ X 



where the path-ordered gauge factors on l.h.s. are omitted for brevity. For the three- 
particle i?-meson DA's the model suggested in [30] is used: 



o 6w o 



\2 

24t^ 



£(7w - 13w + 3£)e 



•2) 



In this model the parameter ljq is equal to the inverse moment of the B meson two- 

(8.3) 



particle DA (j) B . For the latter we use the ansatz suggested in |31 



CO 



1 



2 e °. 9b^) 



. e u 



Appendix C: Functions entering LCSR ( |4.4| ) 

Here we collect the expressions for the coefficient functions Fn D ^ (q 2 ,Lo) entering the answer 



for the correlation function in (4.4) and used to derive the LCSR ( |4.7| ). Below the short- 
hand notation a = u /m B , a = 1 — a is introduced: 



3Q b ( 2m b 

— T {m s - m B a) 

mi \m B 



3Q S 



a 2 m 2 B — m 2 



(2cr - l)(m s - m B cr) , 



JW), 2 \ ?>Qb ( . -, ( 2m b 
F{ '(q ,uj) = — 2-(m s + m B cr) I 1) + 



mi 



m B 



3Q S 



a 2 m 2 n — mi 



: (2a - l)(m s + m B a) 



a 2 m 2 B — m 2 



am B 



Qb 



mt 



2 + 



m b 
m B 



m. 



am B 



3m 2 B Q s 



[a 2 m 2 B — ml) 2 



1 + 


2m b 




m B 


2 




mj 


-a 2 






m% 
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F[ YA \ q \.) 



m? \ m B J a 



4Q S 



2 m 2 B — 'in 



6m s m B Q s 
2 (a 2 m 2 B — m 2 ) 2 



[2a 



\m B J 



Ff A \ q \u) = Fr ) ( q 2 ,u) = 0, 
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III' 



a 



m s \ 2m b ( m s 
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2mb 
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^ ^ m s 4cr — 3 



niB / ms a 



• 2 (2a-l)+m 2 B 



a 2 (2a + 1) + — (2a + (1 - 2a) 



tub 



)— ^) 



^2 w ' w ) = — n 2 i \° 

m( { V m Bj 



+m B (m s - am B ) 
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w-B m B ms' 



+ 



2Q S 



^2~2 2\ m Z \° I 

a z m B — mi. { \ m B J 



m. 



in 



2a 2 + — (4a -1) 



m B 



2q 2 (l - 2a) 



U) 



Finally, we present the substitution relations for the integrals in (4.4): 



o (p — ujv) 2 — m 2 
f2(q 2 ,oj) 



-/ 1 (g 2 , W )e-^ 2 , 



du- 



o 1-u/mB* 



o [(p — ujv) 2 — m 2 } 2 J (\ — bjjm B ) 2 M 2 



+ 



m B (1 - uj /m B ) 2 



f2(q ,0Jq) 



(8.5) 



where f\ : 2 is the product of the function F^ A ^ and the corresponding DA (the latter 
integrated over £), and 



s = umiB + 



m s — q uj/m B 
1 — uj/m B 



1 + 



m 2 — q 2 
(m B - w ) 2 



w 



(so + m 2 B - q 2 ) - yj (s + m| - g 2 ) 2 + 4m 2 B (m 2 - s ) 

2m b 
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